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1 Introduction 

With the start of the operation of the Large Hadron Collider (LHC) in CERN and the 
progress in the analysis of growing data samples for many scattering processes, mastering 
the precise evaluation of strong interactions effects within perturbative Quantum Chro- 
modynamics (QCD) [1-3] will become quickly more and more important. QCD is a well 
established theory - testing the validity of QCD is not an open issue any more. Its principal 
role in the LHC data analysis will be providing precise predictions for rates and distribu- 
tions of quarks, gluons and hopefully other newly found particles carrying colour charge, 
being part of either signal process or background. 

The most important theoretical tools in perturbative QCD (pQCD) calculations, apart 
from Feynman diagrams, renormalization, etc. are the so called factorization theorems, see 
for instance [4-6], which allow to describe any scattering process with a single large mass 
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or transverse momentum scale \ip (enforcing short distance interaction), in terms of the 
on-shell hard process matrix element (ME) squared and convoluted with the ladder parts. 
The hard process is calculated up to a fixed perturbative order. The ladder parts are 
defined for each colored energetic parton entering (exiting) the hard process. The initial 
state ladders are conveniently encapsulated in the inclusive parton distribution functions, 
PDFs. 

The logarithmic dependence of the parton distribution functions (PDFs) on the large 
scale fiF is described as a DGLAP [7] evolution of the PDFs. This evolution was studied 
for the inclusive PDFs up to NLO level in the early 80's, see refs. [8, 9], and was recently 
established at the NNLO level [10, 11]. 

Instead of being encapsulated in the inclusive PDF, the multi-parton emission process 
of the initial state ladder, can be modelled using direct stochastic simulation in terms of 
four-momenta and other quantum numbers, within the Monte Carlo (MC) parton shower. 
Here, the baseline works have been done in mid-80's, see refs. [12, 13], where the leading 
order (LO) ladder was implemented in parton shower MC (PSMC) programs. Standard 
LO level PSMCs implement also the process of hadronization of the light quarks and gluons 
into hadrons and play an important role in the software for all collider experiments because 
of that. 

Fulfilling the challenging requirements on the quality and precision of the pQCD cal- 
culations, needed for the experimental data analysis at the LHC, enforces for the first time 
an urgent solution of the problem of upgrading exclusive PSMC to the complete NLO level, 
the same level, which was reached for inclusive PDFs two decades ago. This task is highly 
nontrivial mainly because the classic factorization theorems [4-6] were never designed for 
the exclusive MC implementation, but rather for defining inclusive PDFs and performing 
fixed order calculations for the hard process, convoluted with these PDFs. 

Let us comment briefly on the longer term physics impact of the present work. Re- 
membering that QCD is not any more a new theory, the main impact of this work will be 
the improvement of calculations of QCD effects, for hadron collider experiments like LHC, 
with the aim of improving chances of discovering directly or indirectly New Physics and/or 
better measurement of the Electroweak Standard Model parameters, especially when high 
statistics, high precision data are accumulated. 

More precisely, this work elaborates on pQCD effects in the initial state, which from 
the perspective of the LHC experiments influence mainly: (A) overall normalization of the 
hard processes through parton luminosities, (B) distributions of transverse momenta {kx) 
of incoming partons deforming many other important distributions in any hard process, 
including searches for supersymmetry, etc., (C) and provide one or more jets accompanying 
hard process. 

The longstanding problem in the data analysis is that the above three classes of 
important QCD effects are addressed by three separate theoretical pQCD calculational 
tools, based on different incompatible perturbative techniques: (A) strictly collinear NLO 
PDFs [7], (B) semi-inclusive schemes of infinite order soft gluon Ut resummation [5, 14, 15] 
or, alternatively, LO parton shower Monte Carlos [12, 16] (C) finite order NLO (NNLO) 
calculations [17, 18] sometimes combined with a LO parton shower MC [19] or collinear 



- 2 - 



PDFs. 

In the analysis of experimental data combining these (and other) techniques is not 
only inconvenient, but also is a serious source of irreducible theoretical (and experimental) 
uncertainties. This old and well known problem becomes more severe with the increasing 
precision of the collider data, and will inevitably plague high statistics, high precision 
LHC data. The ultimate aim of the present work is to provide a basis for designing a 
single Monte Carlo program addressing all three classes of the above pQCD effects at once, 
within the same consistent pQCD theoretical framework. However, for this to be realized 
one has to start with solving one basic difficult problem - the upgrade of the initial state 
parton shower MC to at least the same level as standard PDFs, that is to the NLO level, in 
the fully exclusive manner. The present work provides essential building blocks (exclusive 
kernels) for this critical extension of the parton shower, and analyzes factorization scheme 
differences with respect to the standard CFP MS scheme. While this work focuses on the 
implementation of NLO corrections to the initial state ladder parts (parton showers), the 
hard process part at NLO within the same scheme is discussed in ref. [20] 1 . 

The critical problems to be solved on the way to NLO PSMC are the following: 

1. Violation of four momentum conservation. In the standard collinear factor- 
ization four momentum conservation is broken both between the hard process and 
the ladder as well as between the ladder segments (2PI kernels). The source of this 
non-conservation in the standard collinear factorization is the introduction of the 
projection operators, Pkin m re f- [4], operator F in ref. [9], or operator Z in ref. [24]. 
These operators are absent in the first step of the separation of the collinear singular- 
ities into the ladder parts - they are introduced later on in order to (i) conveniently 
isolate the lightcone variable integration out of the phase space for analytical integra- 
tion and (ii) facilitate the order-by-order pQCD calculations separately for the hard 
process ME and the ladder elements (kernels). This non-conservation happens in the 
transverse momenta, which are anyway treated inclusively (integrated over), hence 
this bad feature of collinear factorization goes almost unnoticed. 2 In the traditional 
LO PSMCs the above non-conservation is repaired "by hand" [12, 13], but in such a 
way that the NLO effects induced by this reparation are analytically almost uncon- 
trollable. This is not a problem, unless one attempts to complete NLO in the ladder, 
or to combine an NLO hard process ME with a LO PSMC [19]. A systematic solution 
of the above problem must involve replacing the projection operator Pkin by a more 
sophisticated operation involving a special parametrization of the entire phase space 
for the hard process and the ladders. An explicit example of such a solution for the 
W/Z production process (Drell-Yan process) can be found in ref. [20]. 

2. Separation of singularities before phase-space integrations. In the collinear 
factorization, separation of the LO singular contributions in the form of the leading 
logarithm In ^ or - poles of dimensional regularization can be done only after the 
phase space integration. In order to construct an efficient Monte Carlo algorithm 

1 A complementary approach can be found in refs. [21-23]. 

2 Except of ref. [24], where it is discussed in a more detail. 
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this separation has to be done at the very beginning, at the integrand level, before 
the phase space integration. This requires going beyond the inherently inclusive 
approach of collinear factorization. The effort of getting the NLO prediction for the 
semi-inclusive distribution in the phase space of the hard process (like kx and rapidity 
distributions in W/Z production) started already quite early, see for instance [17, 25]. 
More recently, Monte Carlo tools combining NLO ME for the hard process with 
the LO PSMC were developed in [19] and [26]. Studies on redefining PDFs in a 
partly exclusive form beyond LO (unintegrated PDFs) [5], or in exclusive form (fully 
unintegrated PDFs) [27] are also pursued. 

3. Negative "probability distributions". The NLO corrections in collinear factor- 
ization are negative in some regions of phase space and therefore cannot be generated 
directly in the Monte Carlo, if we insist on the realistic simulation using positive- 
weight MC events. 3 The reasons for non-positiveness of the NLO corrections is well 
understood. For instance, in the physical gauge positive squares of the Feynman 
diagrams are typically more divergent and form the LO approximation, while non- 
positive interferences are collected in NLO corrections. The use of kinematic pro- 
jection operators in factorization and related subtractions are another sources. The 
factorization procedure collects all these non-positive corrections into separate ob- 
jects and the integration over the phase space is done for each of them separately. 
Consequently, part of the factorization procedure has to be reversed in order to re- 
combine non-positive NLO distributions with the positive LO distributions, before 
the MC algorithm is designed. The above defactorization procedure has to be out- 
lined. An example proposal relying on Bose-Einstein symmetrization was formulated 
in ref. [28] and an even more promissing one is described in ref. [29] (similar to that 
in ref. [30]). 

4. Lack of the published exclusive NLO distributions. We have not found exclu- 
sive NLO distributions forming the NLO corrections in the ladder part in literature 
- all published results in the context of the NLO calculations of kernels for evolution 
of PDFs are integrated over the phase space. The main objective of this paper is to 
provide such distributions for the non-singlet case. 

5. Inclusive treatment of multigluon soft limit. One of the critical issues in any 
type of collinear factorization is the behavior of many-gluon distributions in the soft 
limit. In the inclusive approach it is enough to know that the cancellations between 
real and virtual soft contributions allow us in principle to neglect entire classes of 
diagrams and/or divergent contributions - they sum up to zero. In the exclusive MC 
approach these contributions/singularities, instead of being dropped out, have to be 
modelled precisely to infinite order. The above soft gluon behavior is rather well 
known, albeit more complicated in QCD than in QED due to the (non-abelian) triple 
gluon vertex - the eikonal limit and angular ordering are known to govern it [31]. 
In ref. [32] a detailed analysis of the soft limit for the Cp and CfCa contributions 

3 Unless one admits the painful scenario with negative weight MC events, as in ref. [19]. 
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for non-singlet evolution kernels was performed and the well known angular ordering 
was exposed, both numerically and analytically, for the first time using exact double 
gluon distributions presented explicitly. 

6. Inappropriateness of MS factorization scheme. The issue of factorization 
scheme dependence has to be revisited and one has to decide about the best fac- 
torization scheme (FS) for the PSMC implementation. It should not be taken for 
granted that it will be the MS scheme, which is presently the "industry standard" 
for the inclusive PDFs and in pQCD calculations for the hard process. 4 In the view 
of the above the following question emerges: Can one stay strictly within the MS 
scheme for the MC implementation of the NLO ladders combined with the NLO 
hard process ME? In the present work we will addres this question partly, for the 
MC modelling of the ladder part. It is studied for the hard process ME in ref. [20]. 
Let us indicate that our answer will be negative: FS of the NLO MC has to differ 
from MS scheme for several reasons. Some of them will be discussed in detail in this 
paper. The key element defining FS are the so-called soft collinear counterterms. 5 
In MC we choose them to be identical with the distribution used in the LO MC, e.g. 
the LO level PSMC is constructed by means of iterating soft collinear counterterms. 
This choice on one hand will simplify NLO corrections in the MC, but it will also 
cause departure from the MS FS. The second source of discrepancy will be the fact 
that in the MC, the factorization scale will be identical to a well defined kinematic 
variable Qf of the LO MC, replacing the formal parameter \xp of the MS FS. The 
logarithm of Qf is then used in the MC as an ordered evolution time variable. Qf 
will typically be maximum transverse momentum, maximum angle, or virtuality of 
the emitted particles. For the gluonstrahlung, the choice of the maximum transverse 
momentum results in the MC FS very close to MS FS. 6 However, the soft gluon 
eikonal limit for contributions like gluon pair production or quark-gluon transitions, 
dictates the use of the variable related to the maximum angle (rapidity) as a fac- 
torization scale variable in the MC. Third reason is the presence in the MS FS of 
certain artifacts of dimensional regularization which cannot be implemented in the 
MC in four dimensions. 

7. Constrained evolution. Constrains on the momenta and other quantum numbers 
imposed by the hard process on the initial state parton shower (ladder), especially 
important in the presence of the narrow resonances, must be taken into account by 
the PSMC. The clever and powerfull MC technique of backward evolution of ref. [12] 
is good for the LO PSMC, but for the NLO case one may need something more 
sophisticated. The dependence of the ladder part (PDF) on the factorization scale 
is described in pQCD by the integro-differential DGLAP [7] evolution equation. Its 

4 Since then the FS dependence is reduced to the discussion of the residual dependence on the factorization 
scale [If of MS scheme due to higher orders. 

5 We refer to the exclusive version of soft collinear counterterm, which is the distribution within the 
1-particle phase space encapsulating collinear (and soft) singularity - not its integral as in inclusive FS. 

6 This fact is well known, see also analysis of ref. [33]. 
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solution has the form of a time ordered (T.O.) exponential with the ordering in the 
logarithm of the factorization scale. This T.O. exponential can also be obtained 
directly from the ladder Feynman diagrams. In the Monte Carlo T.O. exponential 
is conveniently modelled using a Markovian MC algorithm. 8 However, a Markovian 
MC algorithm would be highly inefficient for modelling the initial state radiation 
(ISR) ladder, when hard process selects very narrow range of energies and flavors of 
the partons incoming into hard process (with narrow W/Z resonances). Here, the 
backward evolution MC algorithm is a standard solution [12] - used in all standard 
LO PSMCs. It is to be seen whether backward evolution MC is upgradable to NLO 
PSMC. In the meantime, the constrained MC algorithm of refs. [35, 36] offers an 
interesting alternative. While the backward evolution MC needs pretabulated solu- 
tions (PDFs) of the evolution equations, which have to be prepared beforehand using 
non-MC auxiliary codes, the constrained MC performs evolution on its own. 

In this article we mainly addres point 4 in the above list of problems. However also 
point 6 is discussed in some detail. 

Let us now establish precise terminology concerning diagrams and phase space inte- 
gration. We will elaborate on the diagrams contributing to non-singlet evolution kernels 
at the NLO level. 

Generally, in the calculation of the exclusive/inclusive evolution kernels in this work 
we will take paper of Curci Furmanski and Petronzio [9] (CFP) as the starting point and 
as the reference in all NLO calculations. For the non-singlet part of the QCD DGLAP 
evolution we will calculate (or recalculate) both the standard inclusive NLO kernels and 
the new exclusive (unintegrated) ones, which are needed for constructing NLO PSMC. 
This work prepares building blocks for NLO PSMC, whereas the actual MC algorithm and 
all issues related to the factorization scheme used in NLO PSMC will be discussed in a 
separate paper [20]. 

Our main object of interest in the present work are the diagrams depicted in Fig. 1, 
with two emitted on-shell quarks and/or gluons, that is diagrams with two cut lines. We 
will call them 2-real, or shortly 2R, contributions. The other diagrams with 1-real and 1- 
virtual, nicknamed 1V1R, will be also partly considered. Diagrams with 2 virtual (2V) will 
not be discussed, because they will be treated in the same way as in CFP (deduced from 
the baryon number conservation rule, which we keep in the Monte Carlo by construction). 

The NLO 2PI diagrams feature amplitude-squares, depicted in the upper row in Fig. 1: 
double gluon emission 1(a), gluon pair production 1(b) and fermion pair production 1(c) as 
well as interferences, displayed in the lower row in Fig. 1. Interference diagrams enter into 
the MC kernels for the first time at NLO and their correct incorporation in the MC requires 
more care. They can potentially make negative contributions to the kernel, spoiling the 
MC weight. It turns out that the interference diagrams have different statuses in the MC, 
in the following we explain the reasons for it. 

7 Running of the coupling constant is included with the help of the usual renormalization group argument. 
8 Known already since prehistory of the Monte Carlo methods at Los Alamos National Laboratory [34]. 
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Figure 1. Real contributions to NLO non-singlet DGLAP kernel. 



Let us make a simple observation that interferences 1(d) and 1(e) together with squares 
1(a) and 1(b) form the full amplitude square 



r \/\r\ 



+ 



Positiveness of this amplitude square will cause the MC weight introducing these two 
interferences to be positive. Both interferences 1(e) and 1(d) can be implemented in the 
MC in the 2R group of diagrams. 

The interference diagram 1(f) originates from the following amplitude-squared 

where qq pair production 1(c) is already included in the non-singlet class, while qq transi- 
tions amplitude (squared) belongs to the singlet kernel. Apparently, diagram 1(f) corrects 
quark-gluon transitions absent in non-singlet kernel and must be included in the MC to- 
gether with the singlet contributions. 

Similarly, the diagram 1(g) is the interference part in the amplitude square 
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where both amplitudes representing qq transitions essentially belong to the singlet NLO 
kernel. 

In the case of both interferences related to quark-gluon transitions the MC weight is 
not protected by the Schwartz inequality due to the absence of amplitude squares in the 
non-singlet kernels. The above problem will naturally disappear when singlet diagrams are 
included in the game, and we need some temporary fix, while staying in the non-singlet 
class. In order to keep the fermion number conservation we keep these contributions in 
the kernel, but add them to the 1R1V class and treat inclusively. In the following such a 
combination of the 1R1V and integrated subclass of 2R interferences we will call collectively 
unresolved contribution. 

Another important point concerns internal singularities of Feynman diagrams. They 
are present in graphs of Fig. 1(a), 1(b) and 1(c). The double bremsstrahlung diagram 
Fig. 1(a) does not enter into NLO kernel as a whole 9 but only what remains after sub- 
tracting soft collinear counterterm of the CFP factorization. On the other hand, diagram 
of Fig. 1(b) features internal collinear singularity cancelled by the corresponding virtual 
(gluon self-energy) diagram. These diagrams may enter into the unresolved part in the 
MC, or may be modelled in an exclusive manner. In the latter case diagram 1(b) requires 
the construction of a dedicated soft-collinear counterterm which encapsulates the above 
internal singularities and is instrumental in the MC construction. Such a counterterm will 
be defined in Section 3.6. Diagram 1(c) can be also modelled in an exclusive way and also 
needs a dedicated counterterm, see Section 4.3. 

Before that, Section 2 presents an overview of calculations illustrated by the example 
of the LO kernel. Notation and methodology of extracting kernels will be introduced using 
this simple example. The integration procedure for NLO kernels is reviewed in Sections 3 
and 4. The contributions from each appropriate Feynman diagram in the axial gauge are 
calculated in integrated and unintegrated form needed for the MC. Analytical integration 
for control will also be done. Discussion of the collinear and soft singularity structure will 
have the highest priority. 

Sections 3 and 4 summarize results of analytical integrations. In contrast to the results 
presented in [9], 2R contributions will be shown separately instead of the sum of 2R and 
1R1V. Section 5 provides final discussions and conclusions. 



2 Leading order example 



Using the LO diagram of Fig. 2 we are going to introduce some basic notation, which 
will also be useful in the NLO calculations. In addition we will show correspondence 
between elements of the CFP [9] scheme using dimensional regularization (n = 4 + 2e) and 
calculations in the Monte Carlo (n = 4) in this simple example. 

On one hand, the non-singlet bare PDF of CFP [9] , from which the DGLAP kernel is 
extracted, reads as follows 



r LO (x,e) = <5(l-x)(l-Z^) + PP 



[lb 



J dt> n {k) fi- 2e xS(x - ^) C F g 4 W LO (k, e 



(2.1) 



It is not two-particle-irreducible (2PI) 
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Figure 2. Born level diagram. 

where PP is the pole part operator, Cp is the color factor, Zp = 1 + ZM + . . . is the 
fermion wave function renormalization factor. 

Let us explain step by step the other elements in the above formula. One particle 
phase space of emitted on-shell gluon (cut line), in n = 4 + 2e dimensions, reads: 

d*n(k) = -^(2n)5 + (k 2 ) = 1 daa 1+2E rf|a| |a| 1+2 ^ 2+2e . (2.2) 

The emitted gluon 4-momentum k^ is parametrized using Sudakov variables 

k = ap + f3n + k±, (2.3) 

where p is the momentum of the incoming parton and the lightlike vector n defines axial 
gauge. The conditions p 2 = k 2 = lead to relation (3 = — 2a(pn) • Non-abelian coherence 
effects in the soft limit beyond LO are easier to handle if we use the variable 

a= k , (2.4) 

a 

instead of transverse momentum k. Its modulus, |a| = a, will be used to define the 
factorization scale in the MC - we will refer to it as an angular scale. 10 

Phase space in eq. (2.1) requires to be closed from the above, at least temporarily. 
In the CFP scheme this closure plays a marginal role, as e) consists of pure poles. 
The upper limit merely influences intermediate results, through the parametrization of the 
phase space, before taking PP and executing all kind of internal infrared (IR) cancellations. 
On the contrary in the Monte Carlo scheme the variable defining the upper limit of the 
phase space plays an important role of the factorization scale. Its logarithm is the evolution 
time variable in the MC algorithm/code. The most popular choices are: maximum angular 
scale (angular ordering), and maximum transverse momentum of all real emitted partons. 
Virtuality of the emitter line in the ladder or maximum k~ = k° — k 3 are the other valid 
but less attractive choices. 

The function W is just one simple 7-trace factor (see [9]): 

W L °(k,e) = ^~Tr^y^]d a ,(k) = \ 1 + (1 7 2 )2 + £a2 ^. (2.5) 
10 Traditional rapidity variable is equal, up to a constant, to its logarithm, 77 = Mfc Jk ) _ j n j a j _|_ conS £ 
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Setting the upper phase space limit for the angular variable, |a| < Q, we obtain: 

T LO (x, e) = 6 X=1 (1 + zff) + PP{^ (2 ^3+2e / dn 2+2eJ da 8 X - X=a 

2-2a+(l + e)a 2 [°° M |2e _ lfl | 
x al _ 2e ^ ^Ia||a| 2e 0, a |< Q j 

x n i .^h i ^ 2 + 2 - i + ^ 2 - 6(1 - x) 2 g 2 -) (2 - 6) 

= 4=i(i + z F ) + pp| ^ 2e (27r)3+2e (l _ x)1 _ 2e — | 

-A 1 2C F q s ^ 1 + x 2 \ 

~ dx=1+ 2e vr V2(l-^)y + ' 

In the above Zp = — j 2CF K as ( mi — |) provides proper normalization, with regularization 
of the IR pole — > ^z^r^p done exactly as in CFP. The evolution kernel is defined as 
twice the residue of T at e = 0: 

J#W) = 2Res (r^(x, e) ) = (2.7) 

How does the above compare with the Monte Carlo? In the Monte Carlo the same 
integral taken in the limits qo < a < Q in n = 4 looks simpler: 

G LO (x, Q) = (1 - S! 1] )^=i + J d^ A xC F g 4 W LO (k, e = 0)5 x=1 ^ a 6 Q>a>ao 

Q [iu x , 2C F q s / rf 3 fc 1 + x 2 Q 2C F a s ( I + x 2 \ 

= (l-S 1 )5 x=1 + -^- J Wo ^^-^ a >s = ^-—^[ W -^)^ 

Q>a>qo 

(2.8) 



kernel is now the coefficient in front of the collinear logarithm: 



where = In Q 2Cf K as ( In | — |) is the Sudakov formfactor. As we see, the same LO 



»«W> - a^«>(x,«) . ^P(^yj + . (2.9) 
Apparently, 1/e pole of CFP translates into In ^: 

f% a||ari = ^/ Q ^ =ln « (2.10) 

The relation between PDFs and evolution kernels in CFP and Monte Carlo factorization 
schemes is, however, more complicated beyond LO. This is discussed in more detail in 
ref. [37], see also refs. [20, 33]. Generally, all differences between MC and CFP factoriza- 
tion schemes will be traced back to diagrams with subtractions, or with internal collinear 
singularity cancellations. 



1- We will often use shorthand notation 5 a =b = S{a — b), ^2+2 t = r(i+<0 an d Q s = <? 2 /(4tt). 
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3 2R contributions to non-singlet NLO kernels 



In the following we shall calculate bare PDFs and the resulting inclusive evolution kernels 
from 2-real phase space of the Feynman diagrams contributing to the non-singlet NLO 
DGLAP kernels in QCD. Let us stress that our real aim are the distributions over the 2R 
phase space integration. Analytical integrations will be performed mainly as a crosscheck 
with the know results and for testing parts of the Monte Carlo code. 

We shall start with explaining notation and 2R phase space parametrization used in 
the calculations. Differential distributions and 2R phase space integrals will be listed for 
each Feynman diagram separately. 

3.1 Kinematics 

Sudakov parametrization is introduced for both emitted partons: 

h = aip + fan + k i± qi = p - ki fori = 1,2 (3.1) 

where p^ is the momentum of the incoming quark (p 2 = 0) and lightlike n M is the axial 
gauge vector. Real on-shell emitted gluon or quark has 4-momentum kf, and q^ denotes 4- 
momentum of the virtual (off-shell) emitter parton. We will also denote k = k± + & 2 , q = 
p — k. 4-dimensional transverse momenta ki± = (0,kj,0) for i = 1,2 will be also used (to 
be extended to n = 4 + 2e dimensions wherever necessary). From k\ = k 2 = we obtain 

2ai(p ■ n) 

The lightcone variable decreases from 1 to x = 1 — ai — a 2 after two emissions. The angular 
scale variable 

a, = - (3.3) 

cti 

is a preferred choice, instead of transverse momentum. The virtuality of the emitter parton 
after two emissions (entering its propagator) and the gluon pair effective mass are: 

2 2/ \ 2/ \ 1 — a 1 — OL\ o 

q = - aia2 q (a 1 ,a 2 ), g(ai,a 2 ) = a, H a 2 + 2ai • a 2 , 

«2 ol\ (3.4) 

k 2 = «ia 2 a 2 (ai,a 2 ), a 2 (ai,a 2 ) = aj + a 2 - 2ai • a 2 . 

3.2 Inclusive evolution kernels, CFP vs. MC 

In the CFP scheme, the NLO inclusive kernel is extracted from the second order expression 
for the bare PDF, which in compact CFP notation reads: 

^^ .-P^l-'l-PW-^ 1 ^'^ 1 ^^ 1 (3. 5) 

+ Pi4 2] + Pi4 1] ((i - p )4 1] ) + Qpk [ o ] )(pkS- ] ) + o(4) 

where = 1 + Z^p + Zpl is the quark renormalization constant and Kq = + is 
the 2-particle irreducible kernel (truncated to 2-nd order in perturbative expansion) defined 
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M) 



121 

in refs. [9, 38]. The NLO contributions to the evolution kernels are coming from H?Kq and 
Fi^Q 1 ^! — P)ifg ). Diagrams (b-g) in Fig. 1 are in the first class and only diagram (a) is 
in the second class (with subtraction). The other 2-nd order terms like (PKq )(PKq ) and 
Z [ p ] PK [ Q ] yield pure \ poles times elements of LO kernels and do not contribute to NLO 
kernels. The 1-st order lZ [ p ] + P^ 1 was already analyzed in the previous section. From 
now on we drop flavor indices as all but one diagram contributing to the non-singlet kernel 
at NLO level describes qq transitions. The qq flavor indices will be understood implicitly 
if not indicated otherwise. 12 The NLO contribution to the evolution kernel (to bare PDF 
of CFP) from a given 2R Feynman diagram X in Fig. 1 reads: 

P x (x) = 2Res (r x (j;,e)), 

T x (x,e) = PP {^ / d*n(*l)dtt«(*2) xS ( x -p^) Cg 4 W x (h,k 2 ,e) 9 Q>s[ku 

(3.6) 

where the ^-function limits phase space from the above using variable s(k\, k 2 ) = maxjai, 02} 
(resulting CFP kernel is independent of this cut-off), C is a color factor of a diagram X. 

Monte Carlo featuring complete NLO evolution, can be expressed as a time-ordered 
exponential in the logarithm of its factorization scale Q, see [28]. Hence, at the inclusive 
level, it obeys its own evolution equation in InQ with its own inclusive NLO evolution 
kernel, being the derivative in InQ of the MC [28] distribution (truncated to 2-nd order) 13 

d 

d!n~Q 

Here Kq is the same as in CFP and comes from the Feynman diagrams, albeit with real- 
virtual collinear cancellations executed before taking the derivative - so above formula is 
finally executed in n = 4. 

For the use of F' it is enough to apply eqs. (3.9,3.10) below. It acts on the integrand, 
contrary to F of CFP acting on the integrals, hence it provides unintegrated NLO dis- 
tributions for the MC [20, 28]. For our purpose (2R diagrams) the above reduces to the 
following: 



j dLips tfi-^a. V' Q { a K • (1 - P'JM- ( 3 - 7 ) 



Hz) = W ^Q(G b (Q,x) + G a (Q,x)), (3i 



where 



G b (Q,x) = x Jd^ 4 {h)d^A{k 2 ) 5i_ a . =ai+Q2 P'i4 2] 0Q> s (k 1 ,k 2 >qo ( 3 - 9 ) 
contains two-particle-irreducible diagrams and 

G a (Q,x) = x j ^4(^1)^4(^2) 5l-x= ai +a 2 

F'(^ 1] K^)9 Q>s ( kl k2)>qQ - F / (^ 11 )F / (^ 1] )6»Q >s ( fcl)>s(fc2 ) > g 



(3.10) 



contains diagrams requiring soft counterterms. The remaining action of F' is spin projec- 

12 Possible ways of implementing multi-flavor partons in MC have been presented in refs. [39, 40]. 
13 Strictly speaking, diagrams like (b-c) in Fig. 1 produce ln fe (Q//i_R), k > 1 (similarly to higher l/e fc poles 
in CFP), to be resummed into running coupling constant, before applying this formula. 
14 Cut-off go plays no role in evolution kernel. 
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tion, the same way as in CFP. On the other hand, the subtraction term -P' (K^)P' (K^) 
is identical to the double gluonstrahlung LO distribution of the MC, hence it deviates from 
CFP. For example interference diagrams contribute 

G X (Q, x) = J d¥ 4 (fci)d*4(A&) 5^ x=ai+a2 xCg 4 W x (k u k 2 , 0) 9 Q>s(klM)>qo (3.11) 

and subtracted NLO diagrams contribute 

G*(Q,x) = / d9i(ki)d9^(k2) 5i- x=ai+a2 xCg 4 

-i (3-12) 

W x (h, k 2 , 0)0 Q>s{klM)>qo - W ct (k 1 )W ct (k 2 )9 Q>sikl)>sik2)>qo 

Since MC distribution W ct (k\)W ct (k 2 ) encapsulates (by construction) all collinear and soft 
singularities, subtracted G X (Q, x) can be evaluated in n = 4. 

Alternative expressions for NLO inclusive kernels of eq. (3.6) and eq. (3.8) provide 
precisely the same results for graphs (d-e) in Fig. 1, which do not have internal divergences 
nor require subtraction (as in LO case), and are evaluated at n = 4. For diagrams (a-c) 
in Fig. 1 we shall see certain small but important differences, which indicate that the MC 
represents a slightly different factorization scheme than CFP. 

3.3 Overview of the 2R phase space integration 

It is convenient to introduce slightly differently normalized phase space 

j* /-m dU ~ lk 1 dad\-k\.. |2e „ dada. 2 . ; .(2vr) 3 + 2e 

(3.13) 

which is dimensionless in n = 4. 

Most of the presented differential results will be normalized using the above integration 
element. For instance in eq. (3.6) we replace d^ n (k) — > d<& n (k) and 



W 



x 



W*( kl ,k 2 ,e) = C 9 4 x a ^0W x = {^)a\a\a\a\ W X . (3.14) 



(2^) 6 + 4s (2 



Let us outline the general methodology used in the 2R phase space integrations. It 
will be described in n = 4, with some small modifications it will also apply in n = 4 + 2e. 
The integration procedure consists of the following steps: 

(a) Using the identity OQ> ma x{a 1 ,a 2 } = Jo dQ5^ = 

max{ai,a2}' i n t e g ra ti° n variable Q is 

introduced: 

d d Q 

; 1 2 i (3.15) 

da\ da 2 1 f ~ x 

7T / d( P W (01/02,0,0:1,02) O max{a u a 2 }>q o- 

a\ a 2 2tt J 1 AS 
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(b) Dimensionless variables yi = di/Q, Vi 6 [0, 1] are introduced: 



Q _ 

x _ f dai da 2 f dQ f dyi 

J oil a 2 J Q Jo Vl 



'/D 



dm 

Q Jo Vl V2 J 2vr (3.16) 



, Oil, a2)^l=max{j/i,j/ 2 }- 

(c) Integration over overall scale Q is performed: 

Q f dai da 2 f 1 dyi dy : 



x , Q fdaida 2 f dy x dy 2 f 

G = In — / oil — x — ai — a 2 ) I / 

q J ai a 2 J yi y 2 J 

x W x (yi/y 2 ,^ 



2vr (3.17) 



In n = 4 + 2e this integration yields dQQ 4€ 1 — — 



4e ■ 



(d) Nontrivial azimuthal angle dependence enters in the kernels only through the relative 
angle between 2 partons <f> = fa — <p 2 . Integration over <f> and 4> 2 is done. 

(e) Integration over yi, y 2 (eliminating b~i= max { yij y 2 }) is performed. 

(f) Integration over lightcone variables ct\ and. 0,2 (eliminating £=01+02) is 

done using 

the IR regularization of CFP: 15 



L 



1 x da „, . f 1 x da a „. . 
3.4 Gluonstrahlung interference diagram - Bx 

Let us start with the relatively simple ladder interference diagram Bx of Fig. 1(e). The 
expression for 2R dimensionless differential distribution reads 



w B *(k 1 ,k 2 ) = 4c 2 F - 1 -c A c F ) (g) 2 



2„2 



cl-i cL 



l a 2 



q 4 (a!,a 2 ) 



rr^Bx , rrBx & 1 ' & 2 , rpBx & 1 \ & 2 , ^Bx ( & 1 " ^ 

J + J 1 12 hJ 2 ~2 hJ 3 9 2 



a l a 2 a l a 2 



(3.19) 



where: 



Bx _ n * + %2 ( 1 , 1 A ^Bx _ 1 + 2x 2 



T(f x = 2x^^- — + — - 2x, Tf x = — 1 + x-x 

1 — x \ai a 2 J ai 

T 2 Bx = i±^!-i + x -x 2 , Tt = 2(l + x 2 ). 



(3.20) 



Q: 2 

15 This leads to two elementary integrals: I = J* ^ = - In 5, 7i = = — | In 2 <S - , defined 

as in ref. [9]. 
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The first order expression for the PDF in the MC (performing scalar products) reads: 



G Bx {Q,x) 



d* 4 (fel)d*4(fe) a^l-m-oa W'Cfci, k 2 ) 



oo oo 

"l-a—cti+a^ / / i o 

Oil «2 J dl J 0-2 J 27T 





" 2 CaCf J v2^ 



(3.21) 



2„2 



Oi a 



1"2 



T( f - + T B* ^ C os + if* ^ cos </> + if cos 2 , 

Ol\ Cl2 



'max{ai,a2}<Q' 



The above includes factors 1/2! due to Bose-Einstein (BE) symmetrization as well as 2 
multiplying interference diagrams. Following the LO calculation example of Section 3.3, 
the integration over transverse degrees of freedom is done: 



G Bx (Q,x) 



fas\ 2 f da.% da<i 



qo \ 2 / v2jt/ J ai a 2 

2 2 



2x 1 2x(l-a x ) 2 2x(l 



+ if B 



1 



In 



2z(l 
.x 



«2J 



(3.22) 



4~V( 1 -«i)( 1 -«2) 
Integration over a-variables finally provides: 



+ 



aia2 
2x 



G fe (Q,x) = ln — y Bx (x), 



3>~W=(g)-(C|-- Cl <V 



1 



1 + x J 
1 — X 



8I + 8 ln(l - x) -21n 2 (x)l + 4(l + s)ln(x 



(3.23) 



As already said, 7 Bx {x) is the same in CFP and in MC schemes (up to a normalization 
factor 2), 16 because this diagram has no internal collinear divergence. Uncanceled IR 
divergences are still present {Iq term). 

Summarizing, the distribution of eq. (3.19) will enter into the NLO correction to the 
MC exclusive kernel contribution. The distribution of eq. (3.23) will be used for numerical 
overall tests of the MC at the NLO level. 



3.5 Subtracted double bremsstrahlung diagram 

The double bremsstrahlung diagram of Fig. 1(a) (denoted as Br) is not 2PI (it consists of 
two 2PI LO kernels) and needs a subtraction term (referred to as diagram BrC). 

We shall start with a simpler case of integrating Br— BrC contribution to evolution 
kernel in the MC scheme in n = 4. Next we shall recalculate the same Br— BrC contribution 
to the bare PDF and NLO kernel in the CFP scheme, analyzing all its components and 
discussing all differences with the MC case in detail. 



16 There is a difference between normalization in MC and CFP kernels at NLO level 7{x) = 2P(x). It is 
due to the definition of MC kernel as a derivative over In Q not In Q 2 . 
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The differential distributions for two Br diagrams 17 in n = 4 + 2e read as follows: 
W Br (k u k 2 ) = W Brl (ku k 2 ) + W Brl (k 2 , h), 



w B ^(h b \ - 4C f ( as Y a ' a 2 



rpBr I rpl 
1 + J 1 



Br a l - a 2 



+ T Br (e) 



(3.24) 



where: 



rpBr 
J 



T 2 Br (e) = T Br (0) + eTf r , T Br (0) = ^[1 + (1 - a^fx 2 + (1 - ai ) 2 ] 



l + x^ + (l-ai) 2 , Tf 



Br 



-i/_Br rpBr 



,1 - aq 

Ql 
1 



(l + x 2 + (l-ai) 2 ), 



(3.25) 



«7 



[af(x 2 + (1 - qi)^) + a4(l + (1 - ai) 2 )]. 



The most singular term ~ T 2 Br can be rewritten (modulo 0(e 2 ) terms) as a product of two 
LO kernels: 



T Br {e) = ^ (1 _ ai )P^\z 1 ,e)Pi° q \z 2 ,e), 



z,e 



1 + z 2 + e(l - z) 1 
1 - z 



^) + ^ 0) (*), 



(3.26) 



where Z\ = l — a.\ and z 2 = (1 — a± — a 2 )/(l — a±). The above term coincides in the MC 
for the ladder with the following counterterm, being just the LO MC distribution 



K B rO (kl , k2) = ,C\ (0) 2 j^TtiO) = AC 2 ( 



as \ 2 a\a 2 
(1-ai) 2 ^ x ~' ~~ p \2i^) l^ai rq ' ! 



2, 1 _ ai ^ 0> (^° ) (^)- 

(3.27) 

It also enters as a subtraction term into the MC weight which implements NLO corrections. 
The contribution to the inclusive PDF of the NLO MC, including the explicit Bose-Einstein 
(BE) symmetrization factor 1/2!, reads: 

9. <vBr (rr \ _ ! 

q 2! 
rBr/u l. \ _ Tj/Brl/i, 7_ \ , TxrBrli 



G Br (Q,x) = ln^?f: b (x) = ^ I d$ 4 (h)dMk2) S x=1 ^ a2 W B J b 8 Q>max{aha2}>qo , 



WZiki, k 2 ) = W" ri (h, k 2 ) + W* rL (k 2i fci) 

- K BrC (ki, k 2 )6 Q>a2>ai - K BrC (k 2 , h)0 Q>ai>a2 , 



(3.28) 



where 



1 + x 2 



8/ -81n(l-x)+41n 2 (x) 



1 — x 

+ 6(1 - x) - 2(1 - x) ln(x) + (1 + x) ln 2 (x) - 2(1 + x) ln(x 



(3.29) 



is obtained from analytical phase space integration using the same methodology as in the 
CFP scheme, see below. 



Two diagrams because of interchange of vertices due to the Bose-Einstein symmetrization. 
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It should be kept in mind that < ?fub( x ) above is obtained for the angular ordering and it 
would be different if we would have adopted /c T -ordering - the difference would be PKi n ( x ) 
of eq. (3.32) below, see discussion in refs. [33, 37]. 

Let us now turn to the CFP scheme which provides, for this particular diagram, a 
2R contribution to the NLO kernel, different than in the MC scheme. For calculating 
Pi4 1] ((l - V)K [ q ] ) of the bare PDF of eq. (3.5) we follow procedure (a-f) of Section 3.3 
step by step. After integrating over Q, <j) and (f)2, at step (e), we deal in the F(Kq 1 ' Kq) 
part with the singular integral in the variable y = max(yi, 1/2). 

The most singular part due to the \5 y= o term in P(1£q Kq ) is easily integrated: 



1 

+ 2 



1 / Cfcxs\ 



2 



e 2 I n J V- 1 qq ^ 1 qq )\^I\2W 



where (Pqq ®Pqq)(x)\ 2R = ^[41n±+41n(l-x)] + (l+x)lnx-2(l-x) is just the double 
convolution of the LO kernel. The counterterm PKq ((— P)Kq ) contributes the same 
expression, but with (—1) in front, and finally (PK^ )(PK^) adds the same expression 
with (+1) in front. Altogether, the pattern of building correct exponential structure of the 
LO in CFP is much more complicated than in the MC, with a lot of over-subtractions, see 
more discussion in ref. [29]. 

Our main aim however is the residue in front of the - pole. Most of it comes from 

the term (-) in eq. (3.30), which is in close correspondence with the NLO correction 

y ~t 

in the MC. In addition it gets "fall-out" contributions from x e terms. In particular 
rc-independent terms from the expansion 

1 (Q 2 \ 2t ^2e 1 nl ( Q 



+ 2 In + 2oj 2 



(2vr) 2 \fi 2 F J e e \fi 2 F 
luckily cancel between P(Kq Kq ) and the collinear counterterm PKq ((— 1P)Kq ). A sim- 



ilar but x-dependent contribution from T^ 7 " gives rise to a remnant term 
dy 5 -yf J dz 1 dz 2 {^[eP' gq \z 2 )P^\z 1 ) + P^\z 2 )eP^ 



1 

1 f , 5. 



4e 



The second bracket [. . . ] comes from the counterterm - it lacks ^ of the true distribution, 
and one of P'(°)p(°) terms gets killed by the PP operation. Altogether, the above spin 
artifact of CFP (absent in MC) contributes the following: 

Pg r (x) = {^fcl J d Zl dz 2 5 x=ZlZ2 [ - P^) (z2)P (0) (zi) + p(0) (z2)P H0) (zi) ] (3 . 31) 
= (g) 2 ^2(l-x)ln(x). 
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The last important contribution in the class ~ Ji^f^ x e is due to the (aia 2 ) 2e term in 
our particular choice of the phase space parametrization. In fact its role is to cancel the 
dependence on this choice, see also discussion in ref. [33]. It is produced by a similar mech- 
anism of partial cancellation with the counterterm as above, and in the z-parametrization 
reads: 

1 

1 f , & 



4e 



dy 5 -yf j d Zl dz 2 {^[2eHa 1 a 2 )P^(z 2 )P^(z 1 )} - [(2dna 2 )Pff(z 2 )Pj( Zl )]} 

x=ziz 2 

= ±J d Zl dz 2 5 X=Z1Z2 P^(z 2 )P^(z 1 )ln z ^ 1 

Its contribution to the bare PDF and NLO kernel is: 

r£L(z) = ^fln(z), P Kin(x) = (^) 2 C 2 [2(1- *)ln(*)-(l + *)ln 2 (*)] . (3.32) 

Finally the real physics is in the term (^) + in eq. (3.30), which happens to be exactly 
the same (up to the normalization factor 2) as the MC contribution of eq. (3.28). At step 
(f) of the integration procedure it reads: 



r Br _ C F f»s\ 2 f danda-2 f Br aia 2 j 

Lsub - — \2^J J =ai+a2 { ° ~2x~~ 1 



Br a l a 2 



2x(l — OL\ 

ln u 

2 2x{l-a 1 ) 2 ~* ' ±2 2(1 -ai) 2 xa x 



Br al{aia 2 -x) Br a\ (1 - a x ) 2 a 2 

~r J- o — — 7^ — r J o ttt. m 



(3.33) 



After a-integrations we obtain 



Tfub — ^Psub( x )> (3.34) 



where (a) = ^(x) of eq. (3.29). 

Altogether, the CFP kernel from the subtracted Br diagram is: 

P B '\x) = P s B J b (x) + Pgjx) + Pg r (x) 

= ( S) 2 C 2 F f ^ [-4/o - 4 ln(l - x) + 2 In 2 (*)] 



\2ttJ l/ ^T^" rw °"^"^ T m y )i (3-35) 
— — (1 + x) ln 2 (a;) — (1 + x) lnx + 3(1 — x) hxx + 3(1 — x) 

reproducing the result of ref. [9] . 

As noted in ref. [33] Pf[ n (x) is absent in CFP, provided we choose the maximum 
transverse momentum as factorization scale variable: s(ki, k 2 ) = max(A;^, k 2 ). It is simply 
the case because the factor {a\a 2 ) 2e is absent. However, an additional contribution exactly 
equal to P B l n {x) would appear from the 2R integral of eq.(3.28) due to adopting k T - 
ordering, and it would exactly compensate the lack of P B l n {x) from (aia 2 ) 2e . In a sense, 
the CFP scheme features an automatic self- correcting mechanism, such that it provides the 
result for /^-ordering, independently of the kinematic parametrization of the phase space 
actually used. 
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Summarizing on the subtracted Br diagram, the integrand W B £ in eq. (3.28) will enter 
into the NLO correction to the MC distribution, while eqs. (3.29) and (3.35) will be used 
in the numerical tests of the MC codes. The difference between CFP and MC factorization 
schemes for the subtracted Br diagram (for Br+Bx as well) at the inclusive kernel is IR 
finite and reads: 

AcFP-MC^z) = PilnW + P,s p r (x) = 

/rvc\2 (3.36) 
= (|^J C\ [4(1 - x) ln(x) - (1 + x) ln 2 (x)] . 

Let us also sum up C\ contributions of the Br and Bx diagrams (using only C\ part 
of Bx) in the 2R phase space. For the MC scheme we obtain: 

G BrSub+Bx (Q,x) = j d*4(fcl)d*4(*?2) S x=1 ^ a2 (W^ b + W Bx )0Q >laax {al,a2}> qo 

= ln Q rpBrsub+Bx,, 
rpBrsub+Bx^ = yBr (;,.) + yB*( x ) 



V2vr/ 



F 



• } + X ln 2 (x) + (1 + x) (2 ln(x) + ln 2 (x) ) + (1 - x) ( 6 - 2 ln(x) 



(3.37) 



As we see IR part (Iq term) cancels between Br and Bx diagrams. The same phenomenon 
is also true for the differential distributions, see figures and analytical investigation of the 
a.i — > limit in ref. [32], which are based on the above results. The above IR cancellations 
are vital for the stability of the NLO MC weight, see tests of the prototype MC in ref. [28]. 
It should be stressed that it is not guaranteed 18 and here it is true thanks to a good choice 
of the multigluon LO MC distributions compatible with the soft (eikonal) limit already at 
the LO level, see eq. (3.27). 

3.6 Gluon pair production diagram - Vg 

Let us investigate now another important 2R NLO contribution from the gluon pair pro- 
duction diagram Vg of Fig. 1(b). We shall calculate its contribution to inclusive kernels 
focusing on possible differences between MC implementation and CFP scheme. As we shall 
see, the 2R gluon distribution from Vg diagram features very different singularity pattern 
than Br+Bx diagrams discussed previously - it has an internal collinear singularity for the 
mass of produced pair going to zero, a 2 — > 0, that is located in a\ — > 0,2 (instead of a% — > 0). 
Cancellation of this internal singularity happens without an intervention of (1 — F) opera- 
tor, simply by adding a virtual diagram (gluon selfenergy). The remaining residual double 
logarithm (or double pole in e) singularity is related to the running coupling constant. 

Mastering the additional soft-gluon sudakovian IR singularities oti — > will be again 
very important for stability of the NLO MC weight. A complete discussion of the soft 
gluon limit will require introducing the interference diagrams Yg and Bx of Fig. 1(d) and 
Fig. 1(e), hence it will be deferred until the next section. 



3 Such IR cancellations would not work for virtuality ordering in the LO MC, s(ki, k'2) = — (p— k\ — &2 
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From the Feynman diagram we obtain: 
AC A C F 1 / a s V 



T Vg 



2 a \a\ 



1 



T ^ + T V 9{e) ^ + 4 +T V 9 ± 



Vg, 



,2\2 



-2aia 2 + 4(1 - x) - (1 - x){2 - x + x 2 ) I — + — J , T 3 ys (e) = 2x(l + e), 



T^(e) = (l-x) 
T^ 9 = (ai - a 2 ) 



1 



1 



l + x 2 )l^ + 
(1 + x) - (2-x + x 2 ) 



+ 1 
1 - 



(1-x) 



2 2 

a\ a 2 



+ (1 



X 



aia 2 



+ (l + x 2 



(1-x) 

2 2 



(3.38) 



In the above we have kept only those terms 0(e 1 ) from the 7-trace which lead to 

1 1 • (a 2 -a 2 ) 2 

poles, because of the extra - from an internal gluon mass singularity. Using 1 &4 2 = 

1 a2 2 + a4 1 — — — the most singular term in eq. (3.38) is isolated even more 

clearly: 



4CaC f 1 /a s \ 2 afa 2 . 



(1-x) 2 (2tt)^ V2tt 2 / g 4 (ai,a 2 ) 



rnVg . m 
^0 "T J S 



a 2 + a 2 



+ r 



+ lj"(0) 



[(ai + a 2 ) • a] 2 - (a? + a|)a 2 



(3.39) 



The term proportional to T 2 P (e) = T+ 9 (e) + T^ 9 (e) is the only one contributing the double 
pole \ and is explicitly proportional to the product of two LO kernels: 



(1 



+ e 



aia 2 
1 + x 2 «ia 2 
1 — x (1 — x* 



1 + x 2 /a 2 a\ 
1 — x \a\ a<i 



U\Ol2 



1 



X 



+ (1 



where z = a\/{a\ + a 2 ) and P^I, (2) 



a 2 ai 
ai a 2 



(1 



«ia 2 



Pj2W)^J?(z) + 2*e 



(3.40) 



l-z 



+ 



l-z 



+ z(l 



2 • 



What enters into the 2R part of the NLO correction in the MC, see refs. [28, 29], is 
not the above divergent W Vg , but rather the non-divergent difference with the following 
"soft collinear counterterm" (SCC) representing the distribution used in the LO MC: 



K y 3 = (°^\ 2 AC A C F 



2a\a 2 



2„,2 



a, a 



l u 2 



V2vr 2 ; 



aLxa 2 (1 - xY 



r£f(o), 



(3.41) 



where a max = max(ai,a 2 ). It is the result of a slight simplification of the term ~ T 2 _^ 
in eq. (3.39). Generally, such a SCC is not unique, but complete NLO corrections to 
observables are insensitive to its choice. What is highly sensitive, however, is the dispersion 
(and positiveness!) of the MC weight implementing NLO correction. The above choice is 
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unique in this sense, that it encapsulates not only collinear singularity from a 2 — > 0, but 
also all soft gluon singularities Oj — > for all three diagrams Vg+Yg+Bx, see next section. 

P l0] (z) 

This property ensures good behavior of the NLO MC weight. Moreover, the factor 9 ^ 2 
in eq. (3.41) can be iterated into a separate final state LO sub-ladder for the gluon emitted 
from the primary initial state ladder, see refs. [28, 29]. 

As we are working on the exclusive level we are technically similar to the techniques 
of hard process subtractions of ref. [41] (dipole subtraction) or ref. [42] (antenna subtrac- 
tion), 19 but not to the subtractions used in the inclusive calculations of NNLO evolution 
kernels of refs. [10, 11]. 

In view of the above discussion, it is useful to know analytically, for numerical cross- 
check of the MC code and for discussing complete NLO corrections in the ladder MC, the 
following subtracted Vg contribution to the NLO PDF 20 calculated in n = 4 (as usually 
including BE factor) 



G^ b (Q,x) = i| d* 4 (*i)c»4(fc) [w v °-K v ° 



Si- 



x=ai+a2"Q>a ma x>qo 



Q0 



(3-42) 

8/i - 8J + 8 J In ' 




22 

H ln(x) - 81n(l - x) +41n 2 (l - x) - 81n(x)ln(l - x) 

3 3 

Generally, the presence of the SCC subtraction is a natural element in any MC scheme 
with soft gluon (photon) resummation, either to the hard process or to the ladders, in order 
to eliminate possible double counting of the singular term. However, the use of subtractions 
can also simplify non-MC calculations, like analytical integration of the Vg diagram over 
the 2R phase space in n = 4 + 2e dimensions in the CFP scheme, before combining it with 
the gluon self-energy virtual diagram. Let us comment more on that, venturing a little bit 
in the area of combining 2R and 1R1V contributions (complete discussion is beyond the 
scope of this work). In such a case, it is useful to split 2R gluon phase space into a > na max 
and a < Ka max , k <C 1, schematically 

and split the Vg contribution into a subtracted one and the SCC. The subtracted part in 
the decomposition 

v Vg (v Vg _ r CTN , pCT , y CT 

1 V A ) ~r L a>/ta max ' a<fta max 



19 The dipole/antenna subtractions works for at least two ladders (they are dealing with hard process), 
whereas our method works within one ladder. Furthermore, we pay special attention to the fact that our 
counterterms can be iterated in the MC simulation. 

20 In the above l/aj are regularized by a small parameter 5 as in CFP, however, this is not necessary 
once diagrams Yg and Bx are added, see next section. 
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can be evaluated in n = 4, all over the phase space. On the other hand, the SCC part 
T CT is evaluated analytically separately in the "resolved" part a > Ka mSLX in n = 4, and 
separately in the a < Ka max part in n = 4 + 2e. This allows to profit from adjusting phase 
space parametrization to specific complications of the integrand in each part! Finally, one 
combines all three parts into a formula for the 2R integrated Vg in n = 4 + 2e. The 
parameter k, and even the dependence on the particular choice of SCC drops out in the 
final result. The other immediate profit from the above methodology is that one may 
combine 2R from a < Ka max with the 1R1V contribution (gluon self-energy) such that the 
Sudakovian part of the \ pole gets cancelled. 21 

This kind of calculation for Vg in n = 4 + 2e is presented in the Appendix using a 
slightly different choice (for historical reasons) of the SCC: 



Switching from one kind of SCC in the integration to another is relatively simple, see 
Appendix. 

Last but not least, let us discuss the differences between the MC scheme and the CFP 
scheme for the Vg diagram. In the previous case of Br diagram we have seen that the 
basic mechanism of producing differences between two schemes is the action of the 1 — F 
operator. Since this operator is absent for Vg, one generally expects no differences between 
the two schemes. In particular any effect of the terms proportional to e from 7-traces will 
land in the ~ S(a 2 ) part, where soft 2R and 1R1V are combined together in the same way 
in both schemes. 

The only subtle point is the term 60 left over from adding 2R soft and 1R1V con- 
tributions. It comes from integrating over the In — term in the gluon self-energy, and in 
the MC scheme it builds up an a s dependence for some kinematic variable. Which vari- 
able? Changing from one choice of q to another may generate an extra NLO term in the 
kernel (in MC scheme). Our preliminary study shows that taking transverse momentum 
as q is compatible with CFP, that is the Vg diagram contribution is then identical in MC 
and CFP. The contribution from the diagrams Yg and Bx discussed in the next section, 
will contribute the same way in both schemes, due to the lack of any internal collinear 
singularities. 

3.7 Gluon interference diagram - Yg 

The Monte Carlo distribution for the gluon interference diagram of Fig. 1(d) (Yg) is given 




max 




(3.43) 



by: 



W Yg (ki, k 2 ) = W Y9 \k 1 ,k 2 ) + W Y9 \k 2 , h), 




-1 



(3.44) 



21 



The remaining uncancelled part 



C A a s 11 
7T 12 e 2 



is related to the running coupling constant. 
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where: 
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(3.45) 
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2(3 - a?) - 2ai 



2(x 2 -2x + 3) 4(x 3 -x 2 + x-l) 



a i 



The expression for the contribution of the Yg diagram to the NLO kernel is equal to (we 
include both factors 1/2! from BE and 2 due to interference): 



G Yg (Q, 



d^ i (k 1 )d^ 4 (k 2 )W Y 9(k 1 ,k 2 )5 1 . 



x=ai+a2 ^max{oi,a2}<a n 



(3.46) 



The integrand only has the familiar scale singularity which can be extraced in standard 
way as a logarithm. Then the integrals over transverse components take form: 

G^(Q,x) = ln^( 1 -C A C F ) f^-Y / ^^ ( l-x- ai -a 2 ) 
qol — x\2 I \2-n J J a\ a 2 

x / T yg «l«2 _ rjYg Ql"2 T Y g OL\OL2{a\ ~ Ot2) 

\ 2x 1 2x(l-ai) 3 2x(l-x) 



(3.47) 
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After the final integration over on variables we obtain 

Q 



G Y9 {Q,x) = In ^-? Y 9 (x) 
Qo 



\2irJ 



-CaCf 



1 + x 2 
1 — X 



8/i + 16J - 8 J In 



1 -x 



2 In 2 fx) 



- 4 ln 2 (l - x) + 4 ln(x) ln(l - x) - 3 ln(x) + 16 ln(l - x) - 9 - 4Li 2 (x) 



(3.48) 



+ (1 + x)ln(x) + 3(1 - x) + 



1 — x 



Gluon interference diagram Yg features both double and single logarithmic IR diver- 
gences (Ji and Jo). Of course, they must cancel when all diagrams are added. Cancellations 
occur for 2R diagrams not only on the inclusive integrated level but already on the ex- 
clusive unintegrated level, see [32]. We can see them explicitly by adding inclusive kernel 
contributions for gluon interference diagram Yg (3.48), subtracted gluon pair production 
diagram Vg (3.42) and part of gluonstrahlung interference diagram Bx (3.23) proportional 
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to \CaCf colour factor: 



? Y9 (x) + V%(z) + V^ix) = (g) 2 {^CaCf^ 



1 + x 2 
1 - X 



41n(x) ln(l - x) + — ln(x) + ^— 
3 3 



149 



4Li 2 (x) 



(3.49) 



3(1 + x) ln(x) + 3(1 - x) + — + \— — 

1 — x 61 — x 



The above expression for the sum of 2R contribution with colour factor equal ^CaCf is 
free from double and single logarithmic IR divergences (I\ and Jo), as expected. 

4 Other non-singlet diagrams 

The contributions to NLO kernels from diagrams displayed in Fig. 1(c), Fig. 1(f) and 
Fig. 1(g) will be presented in the following. They will be referred to as Vf, Yf and Xf 
diagrams respectively. As discused in the introduction, before all singlet class diagrams 
are included, the contributions from the interference diagrams Xf and Yf should enter the 
MC code in the inclusive form. 



4.1 Interference diagram Xf 

The crossed-ladder diagram Xf contributes the following distribution 
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to the quark-antiquark kernel, where: 
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The integrated distribution is equal to: 
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(4.3) 



and the contribution to the NLO kernel equals: 



G-^(Q,x) = hA^(x). 



(4.4) 



The above agrees with [9] up to the — sign which is a matter of convention. 
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4.2 Fermion interference diagram - Yf 

The differential distribution for the fermion interference diagram Yf of Fig. 1(f) reads: 
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where: 
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(4.6) 



The contribution of the Yf interference diagram to the inclusive kernel (bare PDF) is given 
by (including interference factor 2): 



*"M = (£)>-i<fc°' 
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(4.7) 



This result agrees with that of ref. [43], 22 the difference in sign can be attributed to a 
different definition of space dimension (n = 4 — 2e, e being negative as opposed to CFP). 

4.3 Fermion pair production diagram - Vf 

The fermion pair production diagram Vf of Fig. 1(c) features an internal collinear singu- 
larity when the mass of the produced pair goes to zero. The kinematical structure of the 
Vf graph is quite similar to that of the gluon pair production diagram Vg. The n = 4 + 2e 
dimensional distribution for the Vf diagram reads: 



(l-x) 2 V2 
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(4.8) 



22 The authors use a different regularization technique. In case of the Yf diagram, however, both methods 
give the same result, as explained in [43]. 
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where: 



T Vf (e) = 2(1 + x)ax - 2(1 - x) + 
TP = -2x. 
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The contribution of this diagram integrated in n = 4 + 2e dimensions reads: 
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For exclusive modeling of this diagram we need to deal with its internal collinear singularity 
in a similar way as for the gluon pair production diagram Vg. The following counterterm 
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(4.11) 



can be used both for the MC purpose in n = 4 and for combining real and virtual contribu- 
tions. In the latter case, we decompose the Vf contribution into a part entering Monte Carlo 
(r^ — T v f CT ) + r^/^j^ iax and an unresolved part r^< rea ^ ax required for the cancellation 
of double poles from the virtual contributions. For completeness we give the subtracted 
contribution of the Vf diagram to the NLO PDF: 



d$ 4 (fcl)d$ 4 (fc2) 



5l-x=a 1 +a 2 @Q>a max >qo ~ l n ^ sub^ X ^ 

Q0 



\2irJ 



J 1 



14 
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(4.12) 



5 Summary and outlook 

The main result of this work is a complete collection of 2-real parton (quark, gluon) differen- 
tial (unintegrated) distributions, which enter calculations of the NLO DGLAP non-singlet 
evolution kernels, in a form ready for the use in the Monte Carlo implementation of the lad- 
der (also referred to as NLO parton shower MC). The distributions are given in eqs. (3.19), 
(3.24), (3.38), (3.44), (4.1), (4.5), (4.8). These distributions in the fully differential form 
are not available in the literature. 

We also present the differential distribution of the collinear soft counterterms, which 
are used to subtract internal singularities for some diagrams. These subtractions are also 
present in the MC weights. The MC collinear soft counterterm distributions are defined in 
eqs. (3.27), (3.41). 
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Furthermore, we present analytical integration results. They are presented as the con- 
tributions to evolution kernels from the same 2-real parton differential distributions listed 
above, see for example eqs. (3.37) for all bremsstrahlung diagrams. In case of diagrams 
with internal collinear singularities, subtractions of the MC collinear counterterms is done. 
For certain diagrams it was possible to compare the integration with available published 
results of refs. [9, 43, 44] and agreement was found. 

The QCD evolution of the NLO ladder implemented in the MC is slightly different 
from that of standard MS 1 , as defined and implemented in Curci-Furmanski-Petronzio 
paper [9]. For instace, the differences between MC and CFP schemes are discussed as far 
as it is possible for the 2-real contributions. They are typically present in the diagrams 
with internal, collinear divergences, see for instance eq. (3.36). The complete discussion of 
this issue is beyond the scope of the present work - it will be completed when diagrams 
with 1 real and 1 virtual corrections are added into the game in the forthcoming work. 
Nevertheless, even incomplete results provide us important insight into the differences 
between NLO (integrated) kernels of MC and CFP MS schemes. 23 This analysis will also 
be a practical outcome of the entire project. 

We did not explicitly show results of the numerical cross-checks of the analytical results. 
Let us only mention that all analytical integration results in the paper were cross-checked 
(up to 4-digits) by means of the MC numerical integration using FOAM program [45] within 
the MCdevelop system [46]. 

Summarizing, the present work marks an important step forward on the way to the 
implementation of the complete NLO DGLAP ladder in the Monte Carlo form. 
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A Gluon pair production diagram - Vg 

Let us present details of the calculation of the evolution kernel contribution for the Vg 
diagram. First, we shall show the calculation for Vg subtracted with the counterterm of 
eq. (3.43) in n = 4, then we shall show how to switch to the MC counterterm of eq. (3.41) 
and finally we will integrate Vg over the collinear part of the phase space, a < Ka max , in 
n = 4 + 2e dimensions. 

23 Luckily, all these differences are coming from small subset of diagrams and are relatively simple. 
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A.l Vg - subtracted 

The difference of Vg and the counterterm of eq. (3.43) has no collinear divergence and in 
n = 4 dimensions reads: 
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Finally the ctj integrations result is: 
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A.2 Two counterterms for Vg 

For various purposes it is useful to calculate the contribution from both counterterms of 
eqs. (3.43) and (3.41) in n = 4 with the cutoff |a| > fta max (a > k). Let us start with the 



4 We always implicitly use Principal Value type of regularization for a integrals. 
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Figure 3. The disc represents the region |ai| < |&2 1 - The dimensionless ratio of y = |a_i | / 1&2 1 and 
4>, the relative angle between ai and &2, are shown. Alternative variables (a, 9) are also shown. 



counterterm of eq. (3.43). Introducing dimensionless variables yi = ai/a n 
and performing integration over the overall scale a max we obtain: 



a /flu 



~ a>K ~ %o (1 -x) 2 \2tt) J ai a 2 1_ 



x=a\+ot2 



2vr 



dyi- 



yi 



T 



Vg 



+ 



dy 2 - 



T Vg 



y-2 



q*(l,l)\a 2 (l,y 2 ) 



+ T, 
+ T, 



Vg 



(vl - 1) 2 
■ a 4 (yi,i) 

Vgil-lH 



(A.4) 



a 4 (i,y 2 ) 



'a>K- 



We need to remember that the k cutoff is infinitesimal and all terms 0{n) are neglected. 
Both yi and y 2 integrals are equal because a 2 (yi,y 2 ) is symmetric and q 2 (l,l) = 
does not depend on y^. 
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The integration over a factorizes from the y and <j> integrals, hence it is performed sepa- 
rately. Moreover, because of the cutoff on a 2 it is convenient to change variables. Instead 
of (y, (f>) variables we use (a, 9) depicted in Fig. 3. The jacobian for this transformation is 
equal to 5/vl + a 2 — 2d cos 9. The two integrals above are calculated separately. Firstly, 
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and next 
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Combining both parts together we obtain: 
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After performing the a-integration the final result for the first counterterm of eq. (3.43) 
reads: 
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The integration for the counterterm of eq. (3.41) representing the double gluon distri- 
bution in the Monte Carlo proceeds for the same phase space quite similarly and the final 
results reads: 
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The difference between the two counterterms is, of course, colliner- convergent and it 
reads 
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The above is used to correct the Vg subtracted result of eq. (A. 3) in order to obtain the 
result of eq. (3.42). 

A. 3 2R collinear singularity of Vg in n-dimensions 

For the purpose of combining the 2R contribution with virtual corrections (gluon self- 
energy) it is useful to calculate the contribution to T, bare PDF, or PDF of the MC in 
n = 4 + 2e dimensions, in the collinear region |a| < Ka max . Let us start from the expression 
where dimensionless variables y% = 04/ a max are already introduced and the integration over 
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a-max is performed, giving the - factor: 
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The two integrals over y\ and 2/2 are equal, hence: 
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One more time we see that the a integration factorizes. We start the integration with y 
and <f> and calculate them separately for the 1/d 2 and 1/d 4 parts. We use the variables 
(d, 9) introduced before instead of (y, cf>), then: 
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The e pole is extracted in form of d 2e_1 so we can perform the e expansion of the rest of 
the above expression. We also use an additional approximation, since a is smaller then the 
cutoff k and we are only interested in logs of k we can fix the integration limits of the theta 
integral as and 7r/2, then: 
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7T 



+ e 2d-7rln(2) 



= vr ( - + ln(K) - ln(2)J +0(/c)+0(e), (A.15) 



where we performed expansions in e and k. Performing the same operations for the 1/a 4 
term (the same order for expansions and integrations) we obtain: 



J 2 



= 2 I da f 
Jo Jo 



dy _ 4 9{a < k) 

a 4 {y,l) 

tt/2 



d9a 2e ~ l {2 cos 9 -df 
1 
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27r(-+ln(«)-M2) , 



Vl + d 2 — 2a cos 9 
+ 0{k) + 0{e). 



2 c 



(A.16) 
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Finally, we sum up both contributions and perform the a integration obtaining: 

a<K 2e V2vr/ V 2 

l + x 2 |7 1 . / Q 2 \ , 1\ ( Ar , . 11 



+ ln( ) + 7 -In- )2(4/ + 41n(l-x) 



2e \A-kij, 2 J ' kJ V v 3 

— ln(l — x) 

3 V ; 3 



1 — x 

, s 44, , , 4vr 2 134 

+ 8Ji + 8J ln(l - a?) + 121n 2 (l - x) - — ln(l - x) - — + — 



(A.17) 



/ 11 
+ (l-x) (4/ + 41n(l-x)-y 

The above result will be useful for combining the 2R contribution from Vg with the virtual 
corrections. 
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